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Abstract 

In many applications to biophysics and environmental engineering, sedimentation of 
non-spherical particles for example: ellipsoids, is an important problem. In our work, 
we simulate the dynamics of oblate ellipsoids under gravity. We study the settling 
velocity and the average orientation of the ellipsoids as a function of volume fraction. 
We see that the settling velocity shows a local maximum at the intermmediate 
densities unlike the spheres. The average orientation of the ellipsoids also shows a 
similar local maximum and we observe that this local maximum disappears as the 
Reynolds number is increased. Also, at small volume fractions, we observe that the 
oblate ellipsoids exhibit an orientational clustering effect in alignment with gravity 
accompanied by strong density fluctuations. The vertical and horizontal fluctuations 
of the oblate ellipsoids are small compared to that of the spheres. 
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1 Introduction 



The sedimentation of particles in a fluid under the action of gravity is an 
interesting problem in fluid dynamics and statistical physics which finds ap- 
plication in biology, environmental and industrial engineering. For example, 
chemical reactors, pollution, ink-jet printing, blood, fluidized beds, etc. This 
problem reveals complicated multi-body interactions due to the long-range 
hydrodynamics, that decay for a sphere as 1/r, where r is the distance be- 
tween the particles. The particle velocities strongly fluctuate and there has 
been some controversy about the nature of these velocity fluctuations, [1]. 

In the dilute case and low Reynolds numbers, the behavior of the settling veloc- 
ity V(<f>) as a function of the volume fraction $, in the limit where Brownian 
motion is negligible, has been investigated by Batchelor[2]. V($) decreases 
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monotonically as $ increases following the phenomenological Richardson- Zaki 
law [3]. 



Much of the investigations have been made on the spheres and in a reduced 
manner on slender bodies, e.g. fibers by [1,4]. Fibers have many technological 
applications [7]. An orientational transition is well known for smaller volumer 
fractions, which is characterized by a maximum in the sedimentation velocity 
[8]. Ellipsoidal particles find applications in the models for blood flow [5], car 
paint and sun protectors. 

In our previous work [15], we have constructed a phase-diagram where three 
basic regimes exist namely steady-falling, oscillatory and chaotic regime for 
a single oblate. In this paper, we present results of the simulation of many 
oblates sedimenting in the fluid under gravity, in the steady falling regime of 
a single oblate. The rest of the paper is organized as follows. In section 2, we 
discuss the model, in section 3.1, we present the behavior of sedimentation 
velocity, in 3.2, we discuss the orientation of ellipsoids, in 3.3, we discuss 
orientational order parameter. In section 3.4, we present results at moderate 
Reynolds number with not very large aspect ratios followed by summary and 
conclusion in section 4. 



2 Model 

The model was developed by Hoefler and Schwarzer [10,11], extended by Ku- 
usela et al [12] and has been applied to several studies [14,15]. The fluid mo- 
tion is solved by discretizing on a mesh [11], the incompressible Navier-Stokes 
equations: 



where v is the fluid velocity, p the pressure and / represents an external force, 
which in our problem is gravity. The boundary conditions between the fluid 
and the oblate ellipsoid particles are satisfied considering that the fluid motion 
on the particle surface is subject to the non-slip boundary condition: 



where v t is the translational velocity of the ellipsoid, t(xcm) the vector from 
its center to the point x at the ellipsoid surface, and u> the ellipsoid angular 
velocity. The interaction between the ellipsoidal surface and the fluid adjacent 
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Fig. 1. Snapshot of the oblate ellipsoids falling in a fluid. The picture shows the 
"cluster" formation. The ellipsoid aspect-ratio is Ar = 0.4/1.5, the distance is given 
in units of the larger radius Rm, the Reynolds number is Re = 4 * 10 -2 and the 
particle volume fraction is <5y = 0.05. 



to the ellipsoid is achieved by setting a restoring force that gives rise to a 
"distribution force" in the body term of the Navier Stokes equation. This 
distribution force mimics the presence of the ellipsoids in the sense that the 
fluid inside the ellipsoid moves like a rigid body. A restoring force is applied 
whenever the template particle and the fluid rigid body are not in the same 
position [11]. 

The repulsive force between the ellipsoids is chosen proportional to their over- 
lap. When the oblate ellipsoids are non-overlapping the force is zero and at 
short distances the hydrodynamic forces describing the existence of the fluid 
avoid the contact between the particles [12]. For this force, we choose a con- 
tact function explained in depth in [12] and [13] . The geometry of the oblate 
ellipsoid is characterized by Ar, its aspect-ratio is defined as the ratio of the 
minimum radius R m to the maximum radius, Rm- 

Ar = ^ (3) 

IXm 
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Fig. 2. Mean sedimentation velocity v(t)u for the oblate ellipsoid (dash-squared line) 
and sphere (dash-dot line), as function of the volume fraction <&y. The oblate ellip- 
soid aspect-ratio is Ar = 0.4/1.5, the radius of the equivalent sphere R equ i = 0.97 
and the Reynolds number i?e = 4*10~ 2 . 

We define the equivalence between the sphere and an oblate ellipsoid as the 
sphere that has the same volume, with the equivalent radius: 

Requi — \J RitiRm (4) 



The Reynolds number is defined as: 

Re := 2RMVsPf (5) 



where v s is the mean vertical oblate ellipsoid velocity, 2 * R M being the char- 
acteristic length, pf is the density and v is the dynamical viscosity of the fluid. 
The number of ellipsoids in our simulations is of the order of one thousand. 
We choose the density of the fluid, the Stokes velocity and the larger radius 
of the ellipsoid equal to unity in our system. In all cases the container has a 
height, L = 85 and a base of 12.5 x 12.5, and the lattice constant, h = 0.7. 
The ratio between the density of the oblate ellipsoids and the fluid is 4. 



3 Results 



3.1 Sedimentation velocity for oblate ellipsoids 



In fig. 2 we show the mean vertical sedimentation velocity v(t)\\ as a function 
of the volume fraction $y, within a range of 0.01 to 0.21, for oblate ellipsoids 
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Fig. 3. Velocity fluctuations for ellipsoids (circle-line) and spheres (squared-line) for 
the vertical (solid line), and horizontal (dash-dot line) components corresponding 
to fig. 1. The oblate ellipsoid has an aspect-ratio of Ar = 0.4/1.5, the equivalent 
radius of a sphere is R equ i = 0.97, and the Reynolds number Re = 4 * 10~ 2 . 

and spheres and then compared to the phenomenological Richardson-Zaki law 
= (1— $y) n [3] with n = 5.5. The limit of $y— ► corresponds to the single 
falling ellipsoid which we studied in our previous work [15]. It is interesting to 
point out that the sedimentation velocity of the ellipsoid, is small compared to 
that of the sphere, which follows the phenomenological Richardson-Zaki law. 
This is not the case for fibers (elongated ellipsoids), where it is found that the 
sedimentation velocity has a maximum for smaller volume fraction which can 
exceed the terminal velocity of a single fiber [14]. 

For oblate ellipsoids the mean vertical sedimentation velocity passes through 
a local maximum at $y ~ 0.05. This maximum is quite interesting since 
it is not observed for spheres. Similar non-monotonic sedimentation of non- 
spherical bodies (e. g. fibers) has been reported experimentally by Herzhaft 
et. al.[9] and for prolate ellipsoids in simulations by Kuusela et al. [12] due to 
an orientation parallel to gravity. 

In fig. 3 we present the parallel (||) and perpendicular (_L) components of 
the velocity fluctuations with respect to gravity as a function of the volume 
fraction $y which are denned as: 



The angular brackets indicate averaging over the ellipsoids that have not 
reached the final bottom position at the container. The averages were made 
over at least 50 realizations starting with different random postions and orien- 
tations. In fig. 3 the vertical (parallel to gravity) fluctuations for spheres and 



5V 2 =< Vf > - < V\\ > 2 



(6) 




<V 2 > 



(7) 



5 



16 



-*■ Ratio Sph 
-O- Ratio Ellip 




0.01 



0.06 




0.16 



0.21 



Fig. 4. Ratio of the vertical to the horizontal velocity fluctuations for spheres and 
oblate ellipsoids as function of the volume fraction, $y. The oblate ellipsoid as- 
pect-ratio is Ar = 0.4/1.5, and the Reynolds number Re = 4 * 10~ 2 . 

for ellipsoids are much larger than the respective horizontal components. The 
fluctuations for ellipsoids decrease with the volume fraction. For an equivalent 
system of spheres, the fluctuations show a maximum at intermediate volume 
fractions (<E»y ~ 0.07) [16,17]. In all cases the fluctuations for the spheres are 
considerably larger than the fluctuations for oblate ellipsoids. 

In fig. 4 we present the ratio, SV^/SV^ for spheres and oblate ellipsoids. For 
spheres the ratio shows a maximum around &y ~ 0.07 [16,17]. For ellipsoids 
the ratio shows a slightly larger value than that of the spheres for very small 
volume fractions and has an overall monotonic decrease with the volume frac- 
tion. 

We display the ratio of the vertical velocity fluctuations for spheres to that 
of the ellipsoids in fig. 5. The quotient exhibits a linear behavior with the 
volume fraction following approximately the relation (SV^ ert ^ sph / 5V^ erteUip ) = 
88 * $y + 0.76. The inset shows the horizontal case, also a linear behavior 



In figure 6 we present the parallel, tn and perpendicular, t± components of 
the autocorrelation time for both ellipsoids and spheres. We use the definition 
of the correlation time as: 



Where C(t) is the particle velocity autocorrelation function which is defined 
as C(t) =< SV(t)SV(0) >. Here SV(t) = V(t)- < V > is the local velocity 
fluctuation, where < V > was taken as the mean (horizontal or vertical) 
velocity. 



(SVl rtSph /SV^ elHp ) = 33 * $ y + 3. 



oo 
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Fig. 5. Ratio of the vertical velocity fluctuations for spheres to that of the oblate 
ellipsoids as function of the volume fraction, 3>y. The inset shows the corre- 
sponding ratio for the horizontal fluctuations. The oblate ellipsoid aspect-ratio is 
Ar = 0.4/1.5, and the Reynolds number Re = 4 * 1CT 2 . 
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Fig. 6. Autocorrelation times r (in units of t s ) as function of the volume fraction 
<5y for oblate ellipsoids and the equivalent spheres split into components parallel (| |) 
and perpendicular (J_) to gravity. The oblate ellipsoid aspect-ratio is Ar = 0.4/1.5, 
the equivalent sphere system has a radius R equ i = 0.97 and the Re = 4 * 10~ 2 . The 
results are plotted in a log-log scale. 

The vertical component, shows an overall large value for the oblate ellipsoids 
compared to that of the spheres. The horizontal components of the autocor- 
relation time between oblate ellipsoids and spheres are indistiguishable. 

In figure 6 the four curves decay as power laws given by t pa We see that 
the values of a for the parallel and perpendicular components for ellipsoids are 
cvy ~ 0.11 and a± ~ 0.16 respectively and that for the spheres are ai\\ ~ —0.10 
and a± pa 0.22 respectively 
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Fig. 7. Pair distribution functions for oblate ellipsoids for different volume fractions, 
<3?y. The Reynolds number Re = 4 * 1CT 2 . 

We calculate the pair correlation function for different volume fractions <3V 
and the results are shown in fig. 7. The pair correlation function, for smaller 
volume fractions $y = 0.001, clearly shows large inhomogeneities in the sense 
that there is a "packing formation" as seen in fig. 1, of oblate ellipsoids. These 
inhomogeneities disappear for large volume fraction $y > 0, 2. Furthermore, 
in the intermediate case for <3>y = 0.05 we can see in the pair correlation 
function that the first peak is close to the origin, located at r = 0.6, which is 
also present at $y = 0.2 but smaller. This additional larger peak at $y = 0.05 
could be related to the local maximum in the sedimentation velocity fig. 2. By 
looking at the snapshots as the one shown in fig. 1 one sees that entire bundles 
of aligned particles seem to detach and move down faster which might well be 
the origin of this peak. This kind of "bundle behavior" has also been observed 
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Fig. 8. Euler angles (j), 6 and ip used for the description of the oblate ellipsoid 
orientational behavior [17]. 

in the sedimentation of fibers ref [9] where these bundles settle faster than the 
individual fibers. 



3.2 Orientational behavior 

For the measurement of the orientation we use the Euler angles described in 
fig. 8. The mean vertical orientation (MVO), 9 , as a function of the volume 
fraction, is shown in fig. 9. For smaller volume fraction the MVO shows more 
alignment with gravity and in the limit $y - ► a closer alignment with gravity 
is observed which corresponds to the orientational behavior for one oblate 
ellipsoid observed in [15]. We also see for the MVO an intermediate maximum, 
at $y ~ 0.05, which could explain the local vertical velocity maximum at 
the same volume fraction shown in fig. 2. This intermediate maximum is not 
present for spheres. For larger values of the volume fraction, $y > 0.08 the 
plot shows a monotonic decrease. 

Figure 10 shows the orientational distribution function P(cos(6)) for the ver- 
tical angle, 9, for different volume fractions, $y- For smaller volume frac- 
tions, $y = 0.008 the orientational distribution shows a maximum around 
cos(6) ~ 0.1 in agreement with fig. 9. The limiting case (<E»y — > 0), i.e., one 
sedimenting oblate ellipsoid, studied by us in [15] , presents a vertical alignment 
with gravity (6 ps 90°), and in fig 9 we can see a value of 6 85°. 

As the volume fraction increases, the distributions flatten, and at $y = 0.05 
the distribution shows a moderate maximum around cos(8) ps 0.45, corre- 
sponding to the similar intermediate maximum in figures 9 and 2. We con- 
clude that the vertical velocity is influenced significatively by the orientational 
behavior along gravity, as it is well known for other spheroid systems [8]. 
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Fig. 9. Mean vertical orientation 6 for oblate ellipsoids as a function of the volume 
fraction <3?y. The oblate ellipsoid aspect-ratio is Ar = 0.4/1.5, and the Reynolds 
number Re = 4 * 10~ 2 . 
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Fig. 10. The distribution function P(cos{6)) for the mean vertical orientation 9 
for different volume fractions. The ellipsoid aspect-ratio is Ar = 0.4/1.5, and the 
Reynolds number Re = 4 * 10~ 2 . 

Figure 11 shows the orientational distribution function P(cos(<p)) for the an- 
gle 0, for different volume fractions, $y. The orientation around the vertical 
slightly increases for smaller volume fractions, and decreases with larger vol- 
ume fractions. Similar behavior is also found for the third Euler angle ip. We 
conclude that the Euler angles and ip are not much influenced by the volume 
fraction. 



3.3 Orientational changes 



To quantify the orientation of the oblate ellipsoids we introduce the quan- 
tity \& =< 2cos{9) — 1 > that was also used in [8], [9] as orientational order 
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Fig. 11. The distribution function P(cos(<p)) for the mean vertical orientation <fi 
for oblate ellipsoids. The aspect-ratio is Ar = 0.4/1.5, and the Reynolds number 
Re = 4* 10~ 2 . 

parameter, it would give —1, or +1 if all the oblate ellipsoids were per- 
pendicular to gravity, randomly oriented or aligned with gravity respectively. 
Figure 12 shows the behavior of \& against $y, for smaller volume fractions, 
$y ~ 0.001 — 0.08 the order parameter takes negative values evidencing the 
alignment along gravity and in agreement with the limit, $y — > (one oblate 
ellipsoid) [15]. Approximately at $y ~ 0.08 the order parameter is zero. For 
larger $y > 0.08 a positive order parameter implies the orientation is perpen- 
dicular to gravity. 

In the range of $y pa 0.001 — 0.08, ^ has a local minimum close to <J>y ~ 0.05 
where we found a local maximum in fig. 9 and fig. 2. The simulations were 
repeated with two other different aspects ratios A r = 0.4/0.8; 0.4/2.4 and we 
observed similar behavior. In the case of one oblate ellipsoid ($y — > 0) the 
order parameter \1/ has a value very close to —1 as the ellipsoid aspect-ratio is 
increased [15]. 

3.4 Moderate Reynolds number 

Figure 13 presents the mean vertical sedimentation velocity for oblate ellip- 
soids (□ squared line) and the equivalent spheres (o circle lined) as a function 
of the volume fraction at moderate Reynolds number (Re ~ 7). In our previ- 
ous work this simulation method has been used with success up to Re pa 10 
[11] and [14]. The intermediate maximum for the ellipsoids is not observed in 
fig. 13 as seen in fig. 2 at low Reynolds number. 

A comparison with the phenomenological Richardson- Zaki law (continous line 
in fig 13) shows an exponent around ris p h = 3.2 for spheres and n EUip = 4.0 
for ellipsoids. In both both cases, the data follow the Richardson-Zaki law 
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Fig. 13. Mean sedimentation velocity v(t)\\ for the oblate ellipsoid (dash-squared 
line) and a sphere (dash-circle line), as function of the volume fraction, $y. The 
oblate ellipsoid aspect-ratio is Ar = 0.4/1.5, the equivalent sphere has R equ i = 0.97 
and the Reynolds number Re ~ 7. 

rather closely. These exponents (ns p h = 3.2 and nEiu p = 4.0) are between the 
low particle Reynolds number limit (n ~ 4.5) and a turbulent particle system 
(n « 2.5), [3]. 

Figure 14 (top) presents the vertical distribution function, P(cos(6)) at mod- 
erate Reynolds number. For all volume fractions, P(cos(#)) presents a larger 
distribution around cos(6) ps (6 ps 90°), which tends to be much natter 
(cos(9) > 0.15) than in fig. 10. For the other angular variables, and ip, 
the distributions show a peak around cos(0) w 0, cos(^) ~ 0, and for larger 
volume fractions, they follow a constant behavior. 

The bottom of fig. 14 shows the behavior of the orientational parameter \l/ at 
moderate Reynolds number. For one oblate ellipsoid (<E»y — * 0), the value of 
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Fig. 14. The top picture shows the distribution function, P(cos{8)) of the mean 
vertical orientation at different volume fraction 3>y. The bottom picture shows how 
the order parameter behaves with the volume fraction <&y. The oblate ellipsoid as- 
pect-ratio is Ar = 0.4/1.5, the equivalent sphere has R equ i = 0.97 and the Reynolds 
number Re ~ 7. 

\& is closer to —1 (vertical alignment), as in the case of low Reynolds number 
fig. 11. The intermediate maximum is not seen in figure 12, and the point at 
which the orientational parameter \1/ vanishes, is shifted slightly to the right 
($y ~ 0.1) fig. 14 (bottom). This shift in \1> is also seen in the case of fibers, 
when the Reynolds number increases, by a factor of 5, [14]. 



4 Outlook and Conclusions 



We have simulated the sedimentation of oblate ellipsoids at small volume frac- 
tion ($y < 0.2) and small Reynolds number (Re ~ 10~ 2 ). We have found that 
at intermediate volume fraction the settling velocity exhibits a local maximum 
which to our knowledge has never been reported in the literature. It would be 
desirable to experimentally verify this maximum. 
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This local maximum in the velocity can be related to the non monotonic 
behavior of the vertical orientation of the oblate ellipsoids along gravity, which 
is shown in figures 9, and 10, and can be explained by the "cluster" formation 
shown in fig. 1, which is also found in fiber-like suspensions [9]. 

At low Reynolds number the orientational order parameter \1/ vanishes around 
$y ~ 0.08 fig. 9. As $y decreases the orientational alignment with gravity 
increases as shown in fig. 9 and 14 (bottom), as for low and moderate Reynolds 
number and in the limit $y — > a single ellipsoid aligns with gravity, which 
is a distinctive feature of the steady-state regime for a single oblate ellipsoid 
as reported in references [15] and [19]. 

We also present data at moderate Reynolds number (Re rs 7) for the sedi- 
mentation velocity of oblate ellipsoids as the volume fraction $y is increased. 
As in the case of low Reynolds number the ellipsoids have a smaller sedimen- 
tation velocity than the equivalent spheres, fig. 2 and 13. The data for ellip- 
soids and spheres follow the Richardson- Zaki law [3] with exponents (riEiup ~ 
3.2,Re = 10~ 2 ) and (n Sph w 4.0,i?e = 7) respectively. The P(cos(0)) distribu- 
tion presents a larger alignment of ellipsoids with gravity compared to those 
with small Reynolds number. The vanishing of the order parameter is slightly 
shifted ($y ps 0.1) to the right as the Reynolds number increases (see fig 
14, bottom). The alignment with gravity is present for small and moderate 
Reynolds number as $y — > 0, as shown in fig. 14 (top) and fig. 10, which 
is in agreement with the orientational behavior of a single ellipsoid [15]. All 
the simulations in this work are located in the steady-falling regime of our 
previous work [15] for a single oblate ellipsoid. 

One open problem that should be addressed in the future is the behaviour of 
velocity fluctuations as the container system size and the ellipsoid aspect-ratio 
increases. This work has been carried out for sedimenting spheres but not for 
ellipsoids, [1]. 
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